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Instability of anisotropic cosmological solutions supported by vector fields. 
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Models with vector fields acquiring a non-vanishing vacuum expectation value along one spatial 
direction have been proposed to sustain a prolonged stage of anisotropic accelerated expansion. Such 
models have been used for realizations of early time inflation, with a possible relation to the large 
scale CMB anomalies, or of the late time dark energy. We show that, quite generally, the concrete 
realizations proposed so far are plagued by instabilities (either ghosts, or unstable growth of the 
linearized perturbations) which can be ultimately related to the longitudinal vector polarization 
present in them. Phenomenological results based on these models are therefore unreliable. 



Introduction. Observations of the Cosmic Microwave 
Background (CMB) anisotropics in the WMAP experi- 
ment [l| are in overall agreement with the inflationary 
paradigm. However, certain features of the full sky maps 
seem to be anomalous in the standard picture. These 
anomalies include the low power in the quadrupolc mo- 
ment [3, H, the alignment of the lowest multipoles, 
also known as the 'axis-of-evil' and an asymmetry in 
power between the northern and southern ecliptic hemi- 
spheres [61]. The statistical significance of these effects 
has been debated in the literature. The discussion is com- 
plicated by the difficulty of quantifying the a posteriori 
probability of the effects in the different maps produced 
by foreground cleaning methods and in the context of 
statistical anisotropy where the use of the angular power 
spectrum as a statistic can be misleading. However, re- 
cent studies on properly masked data have shown that an 
anisotropic covariance matrix fits the WMAP \ow-£ data 
at the 3.8cr level Q [l^- The significance of the anoma- 
lous lack of large-angle correlations, together with the 
alignment ofpower has also grown in strength with the 
latest data Q, with only one in 4000 realizations of the 
concordance model in agreement with the observations. 

Such violations of statistical isotropy are considered at 
odds with the standard phase of early inflation. However, 
an albeit more plausible explanation of these anomalies 
arising from a systematic effect or foreground signal af- 
fecting the analysis is not forthcoming. This has led to 
a number of attempts at reconciling some of the anoma- 
lies with the standard inflationary picture through var- 
ious modifications. Specifically, the alignment of lowest 
multipoles could be related to an anisotropic inflationary 
era, whose duration is fine tuned so that the signature 
will be observed in the modes enterin g th e horizon today, 
thus modifying the lowest multipoles [lOl. [Til . [T^ . ITsj . An 
anisotropic expansion has also been considered for late 
time acceleration Although the present statistics of 
the observed Supernovae does not show any evidence for 
the anisotropy [l5|, such studies are motivated by the 
large increase of data that is expected in the next few 
years, with surveys returning many thousands of SNla 
light curves over thousands of square degrees, and by the 
fact that one should keep an open mind on the nature of 
dark energy given our present lack of understanding. 



Anisotropic, but spatially homogeneous backgrounds 
were classified into equivalence classes long ago by 
Bianchi [l^. In the presence of a cosmological con- 
stant and matter fields satisfying strong and dominant 
energy conditions, all Bianchi models, with the possi- 
ble e xcep tion of Bianchi-IX, undergo a rapid isotropiza- 
tion The full formalism for cosmological pertur- 

bations in Bianchi-I backgrounds (the simplest of these 
classes) has been recently developed in [ifl, [H!, [H, HI] , 
with an application to the case in which the only source is 
a slowly rolling inflaton field, which causes the isotropiza- 
tion as an effective cosmological constant. This is the 
simplest anisotropic scenario which can be associated 
with a later inflationary stage. In principle, such a 
background solution can have striking signatures, pro- 
vided that the following inflationary stage is not too long. 
Specifically, different agm coefficients of the CMB multi- 
pole expansion are correlated to each other, and the two 
gravity wave (GW) polarizations behave in a nonstan- 
dard manner and can differ from each other p^ . [ll[ . In 
the case of axisymmetric expansion (equal expansion rate 
in two directions), one of the two GW polarizations expe- 
riences a large growth during the anisotropic era, which 
may result in a large B signal in the CMB [13|. An 
analogous growth is expected also for the general (non 
axisymmetric) case. 

Such simple models, however, do not allow for a small 
and controllable departure from isotropy. Indeed, the 
isotropization due to the (effective) cosmological constant 
starts from a (Kasner-type) singularity [18] and lasts only 
for about one e-fold {At ~ H^^ , where H is the expan- 
sion rate due to the cosmological constant). As a con- 
sequence, one loses predictive power on the initial con- 
ditions for the system. A prolonged anisotropic stage 
can be obtained by introducing some ingredients that vi- 
olate the premises of Wald's theorem ,17J on the rapid 
isotropization of Bianchi universes. This has been real- 
ized through the addition of quadratic curvature invari- 
ants to the gravity action [l^], with the use of the Kalb- 
Ramond axion [23], or of vector fields [2l|. In this Letter, 
we focus on this last possibility, as it is perhaps the sim- 
plest one (at least, from a technical point of view). The 
evidence for an anisotropic covariance matrix reported 
in is based on a primordial power spectrum for the 
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perturbations which is motivated by one of such models 
|22 |; therefore, such constructions deserve close scrutiny. 

In these models, a vector field with non- vanishing spa- 
tial vev is responsible for the anisotropy. To our knowl- 
edge, there are three different realizations of this mech- 
anism. The oldest one dates back to 1989 [2l[, and the 
vev of the vector field is due to a potential V (^^ A^) 
involving only the vector Af^. A more recent proposal is 
characterized by a non- minimal coupling R of the 

vector field to the curvature js^l- For a special value of 
this coupling, the vev of A^ can have a slow roll evolu- 
tion. While the original proposal of this idea [131 realizes 
an inflationary background through several vector fields, 
ref. (25^ suggested a simplified version in which an infla- 
ton scalar field is the main source of expansion, while the 
vector field supports the anisotropy. A completely differ- 
ent class of models makes use of a Lagrange multiplier to 
force a space- like fixed vev for the vector field [22| [ssl ]. 
These three different implementations have been realized 
and studied by several authors [3, [H, [2^, [s^l ■ 

We show that these three class of models contain insta- 
bilities which did not emerge in previous studies. We see 
this from the linearized study of the perturbations around 
the anisotropic inflationary solutions of these models. As 
in all slow roll inflationary backgrounds, each mode of the 
perturbations is initially in the small wavelength regime 
(the wavelength is exponentially small at early times) ; as 
the background inflates, the wavelength becomes larger 
than the Hubble horizon [36| and the mode enters 
the large wavelength regime. This transition is dubbed 
horizon crossing. For the model of [2l|, the system of 
perturbations contains a ghost in the small wavelength 
regime. For the case of the non-minimal coupling with 
curvature of [l^, [IHl , and the fixed-norm case of [22| , the 
ghost appears from some interval of time close to hori- 
zon crossing [37| . The system of linearized perturbations 
blows up at this moment (ssj] . 

Complete computations of cosmological perturbations 
are rather tedious, and the results for the present cases 
can be obtained only through involved algebra. We have 
performed these computations along the lines of [13, llH ■ 
We first write the most general system of perturbations 
(both of the metric and of the vector field); we then 
fix the freedom of general coordinate invariance, we in- 
tegrate out the non-dynamical modes, and finally we 
study the remaining system of dynamical perturbations. 
The divergence of the linearized perturbations is found 
by solving the linearized Einstein equations. Ghost are 
found from studying the kinetic matrix that couples the 
dynamical perturbations in their quadratic action. Such 
computations cannot be reported in this Letter, and due 
to their length, do not provide an insight on the true na- 
ture of the problem. For this reason, we report them in 
a separate and more extended publication [33] ■ The fact 
that the instability is related to the vector field, fortu- 
nately suggests that a partial study, with only the pertur- 
bations of this field included, can shed light on the true 
nature of the problem, without the need to go through 



the technicalities of the full computation. The results of 
this analysis, which are summarized in the next Section, 
show that this is indeed the case. The significance of 
these results is discussed in the concluding Section. 

The instabilities. We assume that the spatial vev of 
the vector field is aligned along the x direction, (A^) 7^ , 
so that the line element is 



ds^ = -dt^ + a {if dx^ + b {tf [dy^ + dz^ 



(1) 



Wc introduce the two expansion rates Ha = a/ a, Hb — 
b/b, and we define their average H and rescaled differ- 
ence h through H = ^"+3^^" and h = Hil-ILl . The 
infiationary expansions that we consider below are char- 
acterized by constant or slowly evolving rates. For the 
models we are considering, h/H = O (-8^) , where B 
is the rescaled vev of the vector field (Ax) = Mp a B 
f2ll m, |2^. Therefore, B must also be slowly rolling 
during the slow roll regime. We consider the phenomeno- 
logically relevant case of moderate anisotropy, B <1 . 

Before studying these models, consider a massive vec- 
tor field in an isotropic background (eq. ([T]) with a ~ b) 
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We assume that has vanishing vev, and we decom- 
pose its fiuctuations as — (ao, diaL + ctf) . The 
transverse vector perturbation af , satisfying diaf — , 
contains two physical modes. These modes are well be- 
haved, and decoupled from the ao, q;l perturbations. 
We disregard them in the following. For 7^ , 
the two perturbations ao, encode one additional de- 
gree of freedom, namely the longitudinal vector polar- 
ization. Indeed the mode ao is non-dynamical, since it 
appears without time derivatives in the action, and must 
be integrated out. Namely, its equation of motion, af- 
ter Fourier decomposition in the spatial directions, gives 
«o I {p^ -l-M^)] ctL , where p = k/a is the physi- 

cal momentum of the mode, k the comoving momentum, 
and dot denotes time differentiation. Inserting this solu- 
tion back into ([2]) we obtain the action for the dynamical 
mode. In Fourier space it reads 
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The longitudinal vector mode exists due to the mass 
term, so it is not a surprise that Af ^ multiplies the kinetic 
term. We see that this mode is a ghost for < . 

Let us now turn to the models of our interest. The two 
models [2l| and [IH can be studied together, using [3§| 
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Expanding the potential at quadratic order in and 
comparing with eq. ([2]) , this action leads to the mass term 
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The equations of motion for the rescaled vev B obtained 
from ^ is 



B + 3HB+QB = 
dV 



(6) 



Q = 2—-2Hh~^h^ - 2h + {I ~ 6^) (2 + + H 

Slow roll of B requires Q ^ (since the 3HB term 
provides a "friction" to the motion). This is achieved 
in two different ways by [2l| and HI 



Ref. |21| studied 
solutions with constant Ha^b in absence of the R term, 
5 = 0. This requires Q = 0, or, in other terms 



dV 



Hh 



2h^ 



-Ha Hb<0. 



(7) 



This corresponds to a negative square mass in cq. ([5]). 
From our discussion of the model ^ we therefore imme- 
diately see that the longitudinal vector polarization is a 
ghost in the limit of isotropic background {B — 0). In 
|25j . the choice ^ = 1/6 is made, so that (following the 
idea of [l^]) the O (H^) contribution is absent from Q . 
Then, slow roll is achieved in the case of small anisotropy, 
S < 1 , and for dV/dA^ < i?^ . We then see that the 
square mass parameter ([5]) is negative in this limit, indi- 
cating that the longitudinal vector polarization is a ghost 
in the isotropic limit in this case too. A more detailed 
study, including also metric perturbations, shows that 
the ghost persists also for moderate anisotropy 13211 . 
Finally, let us discuss the stability of model [23 |. 



S = 



(8) 

In this case, the rescaled vev B is forced to be constant, 
and equal to m/Mp by the Lagrange multiplier A. One 
then finds a background solution with constant expansion 
rates satisfying Hb — (l + B^) Ha We decompose 

the vector field in vev plus fluctuations, ~ (A^) -|- 
(ao, ai, diU + ai) where the index i = 2, 3 spans only 
the coordinates of the y — z plane, and diUi = . The 
equation of motion for A, once expanded at the linearized 
level in the perturbations, reads B ai = Q . This equation 
identically vanishes if the background is isotropic {B = 
0), while it eliminates one of the vector perturbations for 
B ^ . Therefore, contrary to the previous study, we 
cannot consider the isotropic limit in this case. 

We are interested in the quadratic action for the per- 
turbations. It is easy to see that the perturbations a; 
decouple. We are then left with a quadratic action con- 
taining ao and a. ao is non-dynamical in this case too, 
and can be integrated out, leading to the quadratic action 

(5^5 = i y dt d^k a b^p^ {pI -2 Ha Hb) 



(9) 



pi+p^-2HaHb 



where p^ and px are the components of the physical mo- 
mentum along the x— direction and in the perpendicular 



y — z plane, respectively. We see the presence of a ghost 
close to horizon crossing (we recall that Ha and Hb are 
constant, while the physical momentum exponentially de- 
creases, PL oc a~^, pt oc b~^ ). Moreover, the equation of 
motion for a (and the corresponding solution) diverges 
' when the prefactor pj^ — 2 Ha Hb vanishes [40] . These in- 
stabilities are confirmed by the full computation of (32| . 

Discussion. We start by stressing the limits of our 
computation. The above results have been obtained for 
standard kinetic terms for the vector field. Since the U(l) 
symmetry is anyhow broken by the potential term, there 
is no special reason for this choice. Indeed, works on 
lorentz violating vector fields study generalized kinetic 
terms of the type C D V^^'^V^A^ - (32 (V^A^)^ - 
PsW^A'^W i,Af^ . The standard kinetic term corresponds 
to /3i = -Ps = 1/2, ^2 = . Ref. [2l| only discusses the 
case of a standard kinetic term. We have studied pertur- 
bations in this model for arbitrary (3i coefficients. We find 
that the ghost is absent for (32 + 7^ . However, in 
this case one of the perturbations is a tachyon in the early 
time / small wavelength regime [s^ . Ref. [2^ showed that 
the model [2^] is unstable in the case of f3i+ (32 + P3 7^ 
[4l[. For the non-minimal coupling to the curvature, all 
the studies done so far are limited to standard kinetic 
terms, and therefore our computations have also been 
restricted to this case. A second limitation is related to 
the fact that we have performed specific computations 
only for the three models [U [2^. However, these 
models are "prototypes" for the three different ways of 
obtaining the non-vanishing spatial vev explicitly real- 
ized in the literature. We expect that the issues raised 
here also arise in all models for which the anisotropy is 
obtained in one of these three ways. 

Due to these limitations, we do not claim that all possi- 
ble models of anisotropic expansion through vector fields 
are ruled out by our findings. Nonetheless, the issues 
we have found are specific to models with vector fields, 
and should be checked in the stability analysis of all these 
models. Ghosts or tachyons in the early time/small wave- 
length regime are an indication that the vacuum of the 
model is unstable; if these instabilities were present only 
at an energy scale A (inverse wavelength) much greater 
than the bubble scale H, we could simply consider these 
theories as effective ones, valid below that scale A . Cos- 
mological perturbations could therefore be treated as in 
the standard case, and they would lead to predictive re- 
sults that can be trusted up to O {H/ A) corrections (pos- 
sibly, to some higher power). The fact that the instabili- 
ties we have found persist up to horizon crossing indicate 
that A in these cases should be comparable or smaller 
than the Hubble scale. Therefore, even assuming that 
these models can be "cured" , the resulting phenomenol- 
ogy would be profoundly different. Although our con- 
crete studies have produced a negative outcome, we hope 
that some specific constructions will eventually be able 
to overcome the issues found here. 
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